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We propose a simple dynamical model of collective behavior that attempts to 
follow a rule of local neighbor interactions, which was abstracted from exper- 
imental flight data of pigeon flocks. This rule is consistent with the previous 
hypothesis of the basic mechanisms affecting collective motion: short range re- 
pulsion to avoid collisions, longer range attraction to keep the group together, 
and velocity alignment to maintain the same navigational direction. The local 
interactions of our model consist of using naive neighbor estimates: essentially 
assuming that nearest neighbors will move with the same velocity as in the pre- 
vious time interval. The dynamics of our model try to follow the experimental 
rule by using a simple adjustment mechanism with respect to the naive neighbor 
estimates. Prom our simulations we show that by changing the initial conditions 
or the number of individuals involved in the interactions, the model is capable 
of exhibiting a wide range of realistic behaviors, signiflcantly richer than the 
dynamics observed in traditional models. Our study emphasizes the importance 
of using experimental data for making better models of complex systems, and 
that this should contribute to a better understanding of nonlinecir collective 
dynamics. 



I. INTRODUCTION 

Multi-agent dynamical systems with collective behavior have attracted many recent stud- 
ies, especially arising from the scientific interest in animal movement and interactions. Com- 
monly called "swarming", the aggregation and coordination of animals moving together is 
a captivating phenomenon that exhibits intelligent behavior and evolutionary properties, 
which benefit the group as a whole. Such cohesive and synchronized group motion, is essen- 
tial to guide a population through a possibly dangerous environment with many irregularities 
that can definitely affect individual movement. This type of collective behavior has been 
observed and studied in birds, insects, fish, and mammal herds^"^. Research has led to dif- 
ferences regarding the interactions between individuals, in particular, the navigational force 
of a swarm can be dominated by informed "leaders"-, or it could simply be a consequence 
of democratic local interaction rules between neighbors^. The latter is sometimes referred 
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as the "many wrongs" principle, and it asserts that the navigational errors in a group are 
removed by tight cohesion and mutual interaction between individuals^!^. For the particular 
case of bird flocking, recent advances in measurement technology have allowed convincing 
conclusions to be obtained regarding the existence of leaders^, as well as the topological 
nature of the local interactions between neighboring birds^. 

Many types of mathematical models have been proposed to simulate collective dynam- 
ics. Attraction, short-range repulsion to avoid collisions, and velocity alignment, have been 
typically characterized as the main mechanisms to describe swarming behavior in models of 
animal groups^. These rules, influenced by the local interactions between individuals are the 
same principle as that behind the well-known "Boids" model^, which has been traditionally 
used to generate swarm animations. The Vicsek model showed that simple orientation align- 
ment between local neighbors is enough to generate complex collective behavior^S. Newer 
modeling approaches have considered different attraction-repulsion-alignment mechanisms 
and generated a variety of more complex collective dynamics, such as individuals rotating 
around a center like a vortex (also called milling)^"—. Others models have used a set of 
"informed" individuals or leaders with superior information in order to drive the direction 
of the swarnt^ii^'i^. All these models have been based on the principle of manually tuning 
their parameters to get a desired collective behavior in the dynamics, or to compare between 
a wide variety of tendencies in the simulations. New modeling approaches have considered 
positional data to automatically tune parameters in a flxed model structure, in order to 
have behavior based on experimental observations^ii^. The approach we advocate here is 
to begin with actual fleld measurements of animal collective behavior and to ask what rules 
one can deduce directly from that data. The recent availability of extensive experimental 
data motivates this direct analysis, in order to infer the true mechanisms of swarming in 
nature. 

We adhere to the philosophy of using data to construct models capable of simulating 
collective dynamics^, inspired by approaches that have been proposed before in other areas 
and contexts to infer dynamical systems^Si"— , or their natural laws-^. Using a dataset of 
positional data from GPS^, in a previous study we built models that showed the importance 
of mutual local interactions in pigeon flocks^, which is consistent with the "many wrongs" 
principle. From the same study, using our models we inferred an averaged local interaction 
rule that is analogous to the attraction, short-range repulsion, and orientation alignment 
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mechanisms from earlier models. Similar behaviors have been inferred for shoaling fish re- 
cently in another work—. In this paper, we show that by using a new abstracted model in 
which the individuals of the swarm use a simple adjustment mechanism to follow this natural 
rule, it is possible to generate complex and realistic self-organizing dynamics. While many of 
the individual behaviors our model produces have been reported before, such as mobile cohe- 
sive motionii, counter-rotating vortices^""—, initially aligned individuals losing their unified 
movement to uncertain oscillations^^, and pattern formation^^, to our knowledge, it is the 
first one able to generate dynamics resemblant to all of them by simply varying the density 
and the number of interacting neighbors in the simulations. This new variety of behaviors 
emphasizes the emergence of complex phenomena from simple collective interactions and it 
should open the door to a new generation of models with dynamics from behavioral rules 
extracted from actual experimental data. 

II. SEPARATION DYNAMICS FROM EXPERIMENTAL DATA 

In our previous study^^, using a general flocking model built from data of several homing 
pigeon flights'', we measured collective behavior by calculating averaged attraction/repulsion 
statistics that quantify the change in separation of nearest neighbors after a time interval. 
We shall denote this as "separation dynamics" , and describe it by showing how the averaged 
separation between neighbors at a given time, ^{t), affects the change at the next interval 
A^(t). We considered a fixed number of nearest neighbors for interactions, due to a recent 
contribution which showed that birds interact with a specific number of neighbors rather 
than with all within a radius^. From simulations of our models built from experimental 
data, we averaged ^(t) and A^(t) values for different density and speed scenarios^. Our 
analysis produced a curve that showed repulsion tendencies ending at shorter separations 
than 20 m, followed by attraction between 20 and 500 m, reaching a maximum strength near 
120 m. Figure 1(a) shows a cubic spline interpolation of the retrieved curve, which will be 
used in this paper as the main rule for the new simplified model. In addition. Figure 1(b) 
shows speed statistics (magnitude of velocity) also averaged for the same separation values. 
The described behaviors are consistent with the basic demeanors of collective dynamics 
in classical models: repulsion at short separations to avoid collisions, attraction at longer 
separations to keep group cohesion, and velocity alignment to have a synchronized swarm 
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moving together. The latter behavior can be confirmed by visuahzing in the plots that 
near 20 m of separation, we have weak near-zero attraction/repulsion, but a high speed; 
implying steady separations but fast movement, thus equal orientations in the velocities of 
the interacting individuals. We must emphasize that the exact shapes of the curves are 
not important for the new model to be described in the next section, as long as the same 
described qualitative properties (short-range repulsion, attraction, speed) and continuous 
forms are roughly similar. 




^i(t) (m) ^.(t) (m) 

(a) Attraction(-)/repulsion(+) (b) Speed 

FIG. 1: Attraction/repulsion and speed curves obtained from cubic spline interpolations of 
our retrieved statistical^. The plots show change in separation (a) and speed (b) for the 
next time interval {t + At) as a function of average separation to neighbors of an 
individual i at time t, denoted ^(t). Sample time At is 2 seconds. 

III. NAIVE MODEL USING SEPARATION DYNAMICS 

Using the separation dynamics presented in the previous section, we propose to build 
an abstracted 2D model which follows the curves to an extent, especially the attrac- 
tion/repulsion tendencies. Our motivation is to keep the model as simple and tractable 
as possible, in order to present it as a basic reference of how complex behavior can arise 
from relatively simple models that follow tendencies extracted from real experimental data. 
First of all, we shall define the attraction curve as function /[^(t)] and the speed function 
as f [C(^)]; in addition, we denote N individuals, M nearest neighbors, and an infinite space 
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for movement (no boundary conditions). When considering the update for an individual i, 
trying to follow the attraction rule with respect to M neighbors can be quite difficult to 
compute, and has many spatial possibilities. To simplify matters, but still keeping faithful 
to the original concept, we treat the M neighbors interacting with an individual i as a single 
averaged entity, and thus wc shall use the separation of individual i to the average position 
of its neighbors at a time interval, Si{t) — ||xj(i) — {xj{t))M\\, as the separation measure to 
be considered for both curves, and thus we take 5i{t) as ^i{t). We note that even though the 
distance between an individual and the center of mass of its nearest neighbors {Si{t)) is not 
equal to the the average distance between an individual and its nearest neighbors (^j(t)), 
it is a meaningful distance to most of the other neighbors when considering irregularly 
distributed groups, and thus a reasonable simplification. 

The position update for an individual i in the model proceeds in two steps, the first one 
consisting of a preliminary update based on just the speed curve: 

Xi{t + At)' = Xi(t) + v[di{t)]ni{t)At (1) 

where Ui{t) is the unit vector of the previous update, Axj(i) = Xi(t) — Xi(t — At), and 
therefore giving Uj 

(^) = ||Ax-(tj|| - With this first step, wc arc calculating a preliminary 
update which is simply based on keeping the previous orientation (the unit vector), but 
using the new speed value which corresponds to the separation to neighbors at time t, that 
is v[Si{t)]. 

The second step of our position update consists of "adjusting" the preliminary update in 
accordance to the attraction curve. Since we treat the M neighbors interacting with i as a 
single averaged entity, it can be defined as y^it) = {xj{t))M- Taking this into account, we 
define the naive neighbor estimate as: 

y,{t + At)' ^ y,{t) + Ay,{t) (2) 

where Ayj(i) = yi{t) — yi{t — At), which basically is an assumption that all neighbors will 
move with the same velocity as in the previous iteration. We call it "naive" in allusion to 
the traditional "naive predictor" in time series analysis, that consists in using the current 
value as a forecast of the next one. Afterward, we estimate the next separation directly from 
the attraction curve, that is: 

5i{t + At)' = 5iit) + fm)] (3) 
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Then we calculate the "adjustment" , which is a scalar value representing half the distance 
required to reach the estimated next separation, between the preliminary update and the 
naive neighbor estimate: 

<^^it) = + At)' -Mt + At)'|| - S,{t + At)'] (4) 

Finally, we "adjust" the preliminary update and add noise: 

Xi(t + At) = Xi(t + At)' + ai{t)mi{t) + r]{t) (5) 

where mj(t) is a unit vector pointing from the preliminary update to the naive neighbor 
estimate, mj(t) = |jy'(*!|^^j)/Ix'(t+Af)'|| ' ^^"^ Vi^i) is a random vector that represents the noise 
in the update. 

In simple words, the "adjustment" consists of considering all neighbors as a single av- 
eraged entity, and fixing the separation between the preliminary update Xj(t + At)' and 
the naive neighbor estimate yj(t + At)', to attempt to follow the new separation obtained 
from the attraction curve 5i{t + At)'. Of note here is that the | in adjustment aj(t) de- 
notes half of the required magnitude to reach the desired separation, because we loosely 
assume that the averaged neighbor entity will contribute with the other half. The idea 
behind the simplification is to preserve the attraction rule as if individual i is interact- 
ing with a single neighbor (which is actually the center of mass of the neighbors). Under 
steady state conditions, which means M + 1 individuals moving with the same velocity and 
optimal separations (/[5j(t)] = 0), we can expect adjustment aiit) to be zero, since the 
naive neighbor estimate will be accurate, and thus the separation between the preliminary 
update and the naive neighbor estimate will be roughly equal to the desired separation: 
Wy^it + At)' - Xi(t + At)'|| ^ 5,{t + At)'. 

IV. RESULTS 

A. Description of system behaviors 

Extensive simulations using different sets of parameter values were carried out to explore 
the dynamics of this model. We verified that for a fixed number of individuals, the key 
parameters that influence model behavior are the initial positional density (analogous to 
the Visek paper—), the number of nearest neighbors considered in the interactions (M), 
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and the orientations of the initial velocities. For all our simulations, the initial density 
was varied by modifying a parameter Vc, which symbolizes the radius of the circle where 
the initial positions were randomly distributed (uniformly). The noise of the update was 
varied according to parameter e in r]^{t) = U(— e, e). The initial orientations were either set 
randomly or equal (aligned) for all individuals, to compare the differences between initially 
directed and undirected states in the population. 

First of all, some of the most significant behaviors that were observed to be repeatable 
under the same parameter settings, will be described. Afterwards, a general outlook of the 
system behavior for the different parameter values will be outlined as a final picture to 
describe phase transitions and capabilities of the model. All the descriptions and plots are 
restricted to simulations of = 100 individuals, and with At = 2, the latter in accordance 
to the experimental attraction curve. The noise was set to e = 0.5 because it was found to 
be a low enough value that does not dominate the deterministic dynamics of the system, but 
still able to give enough variation for the system to show robust pattern formation and self- 
organization, e.g. individuals "moving around" within a shape. We note that the parameter 
values that we use for the simulations described in this section are not restrictions or strictly 
necessary to reproduce these same behaviors. We chose a parameter space that based on 
our observations of the simulations, showed an interesting range of behaviors and conditions 
that were worthy to analyze. 
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FIG. 2: Simulations with M = 15, Vc = 100, e = 0.5. Random initial orientations in (a), 
snapshot at t = 150. Initially aligned orientations in (b), snapshots at t = 150 and t = 500 

(enhanced online: sim2a.mp4 , sim2b.mp4) 
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For small values of M (between 10 and 20), Tc between 50 and 400, and random initial 
orientations, stable local swarming occurs in small groups, which is expected considering the 
small number of interacting neighbors. Figure 2(a) shows a simulation with M = 15 where 
the population is split into three stable groups that move away. For similar parameter values 
but aligned initial orientations, we consistently observed an interesting behavior where the 
whole population initially self-organizes into a loose triangle, which turns out to be unstable 
and disintegrates as the simulation advances. Figure 2(b) shows an example of how such a 
split can occur. 

Global cohesive movement of the whole population was obtained for values of M > 30, 
but this depended largely on the initial density and orientations. Essentially mid-range high 
density conditions with values near 200, give the best cohesion. For M = 30 and values 
of Tc between 200 and 400, and random initial orientations, the group moves together in 
a single direction after a transient semi-synchronous state, and self-organizes into a stable 
shape that is very resistant to the noise. An example can be seen in Figure 3(a). Using 
initially aligned orientations, high M, and no noise, the population constantly self-organizes 
into more regular shapes with interesting features, e.g. the roughly rectangular, circular, 
and linear sections in Figure 3(b). 
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FIG. 3: Simulations with Tc = 200. M = 30, e = 0.5 and random initial orientations in (a), 
snapshot at t = 1000. M = 50, e = and initially aligned orientations in (b), snapshot at 
t = 1000 (enhanced online: sim3a.mp4, sim3b.mp4) 



A very interesting complex behavior emerges by decreasing the initial density (increasing 
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Tc), and maintaining a high M. Dynamics with particles counter- rotating around a moving 
center, resembhng a vortex, was observed for Tc > 400, M = 50 and random initial orien- 
tations. Figure 4(a) shows a plot of a snapshot during a simulation. The dynamics involve 
individuals further away from the center slowly trying to approach the center with rotations. 
A portion of the population rotates clock-wise around an empty center, while the rest does 
the same but counterclockwise. The empty center of the group, with individuals rotating 
around it, slowly moves in random directions that depend on the overall flow, resembling 
the behavior of many natural collective phenomena. 
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FIG. 4: Simulations with M = 50, = 500, e = 0.5. Random initial orientations in (a) 
snapshot at t = 400. Initially aligned orientations in (b), snapshots at t = 400 and 
t = 1000 (enhanced online: sim4a.mp4 , sim4b.mp4) 



After a long vortex simulation, the group might stabilize into a loosely cohesive unit that 
resembles a comet: a dense core followed by a tail of stranded individuals. This comet-like 
behavior can be obtained easier in a different simulation, by using similar parameters but 
considering initially aligned orientations instead. Figure 4(b) shows such an illustration of 
how the dense core is formed early on, and how the following individuals form a "line" which 
lags behind, resembling a comet and its dissipating tail. 

Under the same conditions as the comet dynamics (high M and aligned orientations), 
decreasing the initial density produces an interesting behavioral transition. Instead of going 
from uncertainty to order as in the vortex simulation (initially oscillating movement that 
later stabilizes), initially aligned individuals will later on reach opposing directions in the 
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core as they try to attract each other, and form a counter-rotating center reminiscent of the 
vortex previously described: see Figure 5. Once the second state is reached, the population 
loses its initial global direction, and is anchored around the vortex center which can slightly 
move around randomly. The overall behavior can be described by saying that what would 
have been the core of the "comet" becomes a vortex instead, causing oscillations. From our 
observations, the system consistently has this initial transition from order to uncertainty 
under these conditions, and it usually does not stabilize into global directed movement. We 
will denote this state as the semi-vortex. 
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FIG. 5: From order to uncertainty (semi- vortex): a simulation with M = 50, Vc = 1000, 
e = 0.5, and initially aligned orientations, at t = 500 (enhanced online: sim5.mp4) 



Changing to high density simulations, we observed interesting "convergence" behaviors 
for some parameter ranges and alignments. Basically, two different groups are formed ini- 
tially, resulting from the short-range repulsion at high densities, but later on they converge 
smoothly and re-organize into a single unit. This behavior was observed very frequently 
for M = 50 at very high densities [r^ between 10 and 50), and random initial orientations; 
though it was more consistent if aligning each half of the population with a 71/2 angle differ- 
ence. Figure 6(a) shows a snapshot of a simulation that exhibited this convergence, which 
results from a weak interaction still influencing each group due to the large value of M, and 
causing a future re-grouping. 

Finally, for the similar parameters but with aligned orientations, the individuals form a 
compact cohesive unit and move together, as can be observed in Figure 6(b). Of interest 
here is that the individuals move within the pattern without affecting the shape or the global 
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FIG. 6: Simulations with M = 50, Vc = 50, e = 0.5. Divided initial orientations in (a), half 
the population pointing at and the other half at tt/2, snapshots at t = 500 and t = 1200. 
Initially aligned orientations in (b), snapshot at t = 500 (enhanced online: sim6a.mp4 , 

sim6b.mp4) 



direction of the group, closely resembling the behavior of a bee swarm. This is likely due to 
the short-range repulsion and noise, creating sensitive deviations to a congestive formation. 
This behavior is consistent with the "many wrongs" principle-, in the sense that the many 
errors, deviations, or noise, within the individuals, are averaged out to form a single cohesive 
global direction which influences the navigation of the whole group. 

In order to give a better general picture of how changing density and speed affect the 
dynamics, we performed extensive manual visualization on simulations and classified them 
according to some of the behaviors just described. For this task, the initial radius Vc was 
varied from 100 to 1000 in intervals of 20, and M was spanned from 5 to 50 in unit intervals; 
thus giving 46x46 different parameter scenarios. These particular ranges were chosen because 
they were enough to span the significant behaviors in the model. To decisively associate a 
particular parameter setting to a behavior, we decided that at least four out of five of its 
simulations must show the behavior. The cases with two possible behaviors with similar 
probabilities, were classified as "variable behavior" . The results from our visualizations can 
be seen in Figure 7, for both random and aligned initial orientations. We must emphasize 
that due to the manual visualization in the extensive number of simulations considered 
(large number of parameter settings, and the repetitions to ensure proper classification), 
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the behavior zones are only approximate. Nevertheless, they give a good picture of how 
certain behaviors can be obtained and the phase transitions by changing either density or 
the number of interacting neighbors. 
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FIG. 7: Rough sketches of behavior zones for different initial densities (radius Tc) and 
interacting neighbors (M); with a resolution of 46x46 parameter settings. All simulations 

had N = 100 and e = 0.5 



Figure 7(a) shows that for initially random orientations, a transition from comet-like 
behavior to vortex can occur by increasing M at low densities (r^ > 450). The same 
transition occurs by fixing M near 40 and increasing Tc at low densities. Cohesive behavior 
is associated with higher densities, and its change to either comet-like behavior or a vortex, 
depends of the value of M when decreasing density; again a larger M will cause vortex- 
like behavior. When considering initially aligned orientations in Figure 7(b), we have the 
change between unstable global movement to a cohesive unit by either increasing M at high 
densities, or decreasing the density for values of M near 30. The latter implies that higher 
densities do not necessarily imply best cohesion. The semi-vortex behavior only occupies 
the far corner region of low density and high M, which tells that it will only emerge in 
cases of high separation and high degree of interaction. As a final but very interesting 
note, it is surprising to see that for near 350 and large M, initially random orientations 
give cohesive behavior, but aligned orientations produce comet-like behavior with stranded 
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individuals. This occurs in the latter because the initial equal orientation in the group 
causes the individuals at the front to synchronize fast as a single unit, while the rear part 
of the population struggles to keep up due to the weak attraction and low speeds for longer 
separations (see Figure 1 curves). This is actually the essence of the comet-like behavior, 
since the weak attraction is what causes the "tail", as can be seen in figure 4(b). When using 
random initial orientations, the longer initial synchronization time actually helps the group 
stay together, and self-organize into a stable structure. It may be that the shapes formed 
by an initially disordered state are more stable than the ones formed with equal initial 
velocities, which underlines the complexity of the model dynamics and their structure. 



B. Quantitative measures of order 



Finding numerical measures that can quantify all the behaviors we have presented is a 
difficult task, and this is why we presented a phase diagram based on manual visualization. 
Nevertheless, we now use two order parameters to identify the global alignment and "vor- 
ticity" in the system. First of all, the magnitude of the average velocity of the system at a 
given time interval: 



1 ^ 

i=l 



(6) 



is a good measure of the global order or direction of the system. A disordered system with 
individuals moving randomly should give a near zero Vm{t) because the velocities cancel out, 
while an aligned swarm results in Vm(t) equal to the average global speed of the group. The 
same or very similar measures have been used in previous studies for the same purpose^i^i. 
To quantify the "vorticity" of the swarm, i.e. how close it is to a vortex state, we measure 
the average absolute value of the angular velocities of the swarm at a given time: 

1 ^ 

wm(t) = ^5^|c^^(t)| (7) 

j=i 



which tells us how much the individuals are rotating around the center of mass of the swarm. 
We chose the average absolute value {j^J^iLi l'^j(^)l) instead of the absolute value of the 



average 



k Sj=i ^ii^) ) '^^ angular velocities, in order to be able to identify the counter- 
rotating vortices, since the angular movement of many individuals would be opposite and 
cancel out with the latter measure, giving theoretically a low value which would make it 
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FIG. 8: Averaged Vm{t) and ujm{t) values for 50 simulations of 1000 time intervals at 
different initial densities (radius Vc) and M = 50. Plots (a) and (c) correspond to random 
initial orientations, while (b) and (d) to aligned initial orientations. The transitions from a 

cohesive state to a vortex can be seen by the decrease of global direction Vm{t), and 
increase of oOmit)- The vertical lines correspond to the manual observations from Figure 7. 

difficult to distinguish from a cohesive state. In general, we expect Umit) to be close to 
zero for an aligned noise-less flock heading in the same direction, but obviously greater for 
a system in a vortex state (regardless of the direction), with particles rotating around. The 
angular velocities uji{t) are estimated by taking phase differences: uji{t) = 9i(t) — 9i(t — 1), 
where di{t) is calculated from the x-axis to the relative position of individual i to the mean 
position of the group at time t: Xi(t) — (x(t)). 

With the purpose of verifying numerically the transition from a cohesive swarm to a 
vortex, we fix M = 50 and span the same region of space as considered in the previous 
subsection, since we can see in Figure 7 that M = 50 is a good enough value to cover the 
transition between these behaviors. For each Tc value from 100 to 1000 (with increments 
of 20), we averaged Vm(t) and 0Jm(t) over 50 simulations and 1000 time intervals. Figure 
8(a) shows that for random initial alignments, Vm captures the loss of global direction as 
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we go from the cohesive to the vortex state at around rc = 400. The increase in vorticity 
can be seen in Figure 8(c) for Um, but interestingly, we see a drop near Tc = 1000, and 
this imphes that at lower densities the vortex loses force because some individuals do not 
interact as much. For initially aligned individuals. Figure 8(b) shows how the decrease 
in global direction is less strong, but decreases more significantly when we approach the 
transition from the comet to the semi-vortex state. This tells us that as we leave the comet 
state, the swarm loses global direction and the uncertainty in the system increases. Again 
we confirm this in Figure 8(d), with a very interesting reaction of Um as we approach the 
buffer area between the comet and the semi- vortex. The vorticity measure slightly increases 
but becomes very noisy as we approach = 700, and this implies that the global behavior 
is very uncertain in this parameter region. 
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FIG. 9: Time courses of Vm.{t) and ujm{t) for three scenarios (M = 50): cohesive unit 
(random initial orientations, r^. = 200), vortex (random initial orientations, = 800) and 
semi- vortex (aligned initial orientations, Vc = 950). The time series were averaged over 50 

simulations each. 



In order to characterize the time evolution of the phase transitions. Figure 9 shows the 
averaged time courses of Vmit) and Umit) for simulations that result in cohesive, vortex and 
semi- vortex states. In Figure 9(a) we can see how both the cohesive and vortex simulations 
start at disordered states (low Vc(t)), with the former steadily self-organizing into a cohe- 
sive flock and a stable Vc{t), and the latter reaching no consensus. The semi- vortex runs 
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start with an initially aligned group (high Vc{t)) which loses its direction as time advances. 
Figure 9(b) shows how the cohesive simulations start with a very high vorticity rate in the 
synchronization stage, and later it is decreased as the swarm aligns. Both the vortex and 
semi- vortex simulations have increases of Um{t) as time advances, with higher values for the 
former. Interestingly, the measures tell us how the cohesive simulations start in a disordered 
state and later self-organize into a stable group, as well as how for semi-vortex simulations, 
an initially aligned swarm evolves into an uncertain state that loses its global direction and 
results in some individuals forming a smaller vortex. In general, the transitions in Figure 
9 are consistent with our observations of the previous subsection. These final illustrations 
show us how both order parameters provide a clear description of the global alignment and 
vorticity of the population, and can characterize the transitions that occur. 

V. CONCLUSIONS 

We must emphasize that roughly the same behaviors as the ones described were also 
observed qualitatively for swarms of smaller or larger sizes. As should be expected for 
different population sizes, the parameters required to get a desired behavior are different, 
due to the difference in density properties and neighbor interactions. Nevertheless, even 
though we tried to summarize the most significant behaviors that were found in our extensive 
simulations, there might be other relevant behaviors that we have not discovered and require 
more exhaustive computation or a deeper mathematical analysis. Still, we present our results 
as an initial taste of the wide variety of realistic behaviors that can be obtained from our 
simple model based on attraction and velocity distributions abstracted from pigeon flock 
data. 

In our previous work^^, we built a computational model from experimental data and 
demonstrated that it is capable of performing accurate simulations of collective pigeon flights. 
In this paper we used a simple attraction/repulsion rule extracted from this model. We find 
that this rule is similar though significantly different from those rules used in the standard 
models of collective behavior-*^. Our main motivation in this current work is to make a 
simpler abstracted naive model: a model based on this inferred rule, which is nonetheless 
tractable and capable of simulating a wide range of phenomena. We then simulated this naive 
model on a range of parameter values far broader than the original experimental system, 
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and demonstrated it is capable of emulating a considerable variety of realistic collective 
behaviors that resemble natural phenomena. 

Several models have individually reported similar behaviors as the ones we have described. 
To name a few, Couzin et al. documented swarming, cohesive, and milling (vortex) dynam- 
ics by varying the interaction zones of individuals^^. Chuang et alM and Touma et al.— 
have shown counter-rotating vortices by changing the attraction and repulsion strengths of 
their models. Similar dynamics to our semi-vortex state, where initially aligned individ- 
uals later produce oscillations, have been described by Erdmann et al.— , using a model 
which only considers attraction between individuals and noise. The strength of our model 
lies in the fact that by following a rule of the basic attraction and repulsion mechanisms 
extracted from experimental data, we can generate this wide range of interesting complex 
behaviors. In addition, to our knowledge, the comet dynamics (Figure 4(b)) have not been 
documented before. This behavior occurs because the weak attractions at longer separa- 
tions in our curve (see Figure 1(a)) cause the more distant individuals to lag behind the 
core in low density scenarios, and thus are a result of using our experimentally motivated 
attraction/repulsion dynamics. Pattern transitions are also now being considered in models 
for collective dynamics^, and furthermore, our simulations also show pattern formations 
such as loose triangles (Figure 2(b)) and other shapes (Figure 3(b)), making the range of 
behaviors even larger. Finally, we stress that our model is capable of producing complex 
self-organized behavior without a global navigational mechanism such as a "leader" , exter- 
nal perturbations from the terrain, or a particularly fixed particle distribution; and this 
is achieved by only changing the initial density conditions and the number of interacting 
neighbors. 

Also worth mentioning is the fact that the stable behaviors observed in our model show 
higher densities at the shape borders than in the center. This type of density distribution 
is consistent with what was measured in a recent investigation of starling flocks from real 
dataM: The existence of optimal density conditions, implying that initially high densities do 
not guarantee best cohesion, which was observed in our previous study^, was reconfirmed 
with our model. Our simulations are also consistent with the "many wrongs" principle^, since 
in many of the runs that showed cohesive behavior, the individual errors of the individuals 
did not affect the global navigational direction of the whole swarm, and it showed strong, 
cohesive, and stable motion. We believe that what is presented in this paper is simply the tip 
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of the iceberg of what could be obtained by building mathematical models which attempt 
to follow or preserve rules extracted from real experimental data, especially in collective 
systems. 
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